RANDOM POINT SETS AND THEIR DIFFRACTION 

MICHAEL BAAKE AND HOLGER KOSTERS 

Abstract. The diffraction of various random subsets of the integer lattice Z d , such as the 
coin tossing and related systems, are well understood. Here, we go one important step 
beyond and consider random point sets in R d . We present several systems with an effective 

CN| , stochastic interaction that still allow for explicit calculations of the autocorrelation and the 

diffraction measure. We concentrate on one-dimensional examples for illustrative purposes, 
and briefly indicate possible generalisations to higher dimensions. 

In particular, we discuss the stationary Poisson process in R d and the renewal process on 
the line. The latter permits a unified approach to a rather large class of one-dimensional 
structures, including random tilings. Moreover, we present some stationary point processes 
that are derived from the classical random matrix ensembles as introduced in the pioneer- 
ing work of Dyson and Ginibre. Their re-consideration from the diffraction point of view 
1 1 improves the intuition on systems with randomness and mixed spectra. 

■*— > 



3 

On 



1. Introduction 

Mathematical diffraction theory is an abstraction of kinematic diffraction |10j that is both 

mathematically rich and practically relevant. Its insight helps to explore the true setting 

and difficulty of the inverse problem of structure determination from (kinematic) diffraction 

data; compare [3[ [5j [16] and references therein. Moreover, mathematical diffraction theory 

£»-. ■ has several important connections with the theory of dynamical systems [U El IH [9] . 

While this is well studied in the classic case of crystals and, more generally, pure point 

diffr active systems, much less is known about structures with continuous (or mixed) diffraction 

spectra. The Thue-Morse chain is a well-known example with singular continuous diffraction, 

but mixed dynamical spectrum; see [H \TT\ [26] for more. Recently, some progress has also 

^ ■ been made for point sets of stochastic origin; see [2] and references therein for a survey. 

Nevertheless, the collection of fully worked out and understood examples is relatively mea- 
gre in comparison with the situation of pure point diffractive systems. Here, we summarise 
two important and versatile examples from [2j and augment them with two examples from the 
theory of random matrices. The latter essentially derive from old papers by Dyson [12| and 
Ginibre [15j . which were later re-analysed by Mehta [21], though the results seem unknown 
in diffraction theory. 

Let us briefly recall the setting of mathematical diffraction theory; see |18[ [2] and references 
therein for more. The underlying structure is modelled by an essentially translation bounded 
measure u, which may be signed or even complex; see [22] for background on measure theory. 
The corresponding autocorrelation measure 7, or autocorrelation for short, is defined as a 
volume- weighted limit, 

1 7 = lim -tttt-t , 
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where uj r denotes the restriction of oj to the (open) ball of radius R around the origin. More- 
over, ui is the 'hipped-over' measure defined by uj(g) = uQj) for arbitrary continuous functions 
of compact support, with g(x) = g(—x). We implicitly use the Riesz-Markov representation 
theorem that allows us to identify regular Borel measures with linear functionals on the space 
of continuous functions with compact support. In general, the limit in ([T|) need not exist, 
but we will only consider situations where it does, at least almost surely in the probabilistic 
sense. 

By construction, 7 is a positive definite measure, hence it is always Fourier transformable. 
Here, we follow the convention of [2J, with the factor 2tt in the exponent (via e _27nfcx ) rather 
than in front of the integral. The result is the diffraction 7, which is a positive measure (by 
the Bochner-Schwartz theorem [22]) that describes, loosely speaking, how much intensity is 
scattered into any given volume of space; see [X0|. [18] for more. The diffraction measure has 
a unique decomposition 

( 2 ) 7=(7) pp + (7) sc +(7) ac 

into its pure point part (which is a countable sum of Dirac measures, known as Bragg peaks), 
its absolutely continuous part (which comprises everything that can be expressed by a locally 
integrable density relative to Lebesgue measure, known as diffuse scattering) and its singular 
continuous part (which is everything that remains, and is often disregarded in crystallography 
- but see [2H]h 

The focus of this paper is on systems with structural disorder, which means that we will see 
either purely absolutely continuous spectra or mixtures thereof with pure point spectra. We 
begin with a review of the Poisson process and the renewal process, where we follow [2] and 
adapt it to the concrete setting of crystallography. Then, we consider some point processes 
that can be extracted from random matrix theory and give rise to further examples that are 
explicitly computable. 

2. Poisson process 

The homogeneous Poisson process in R is an ergodic point process that is often considered 
as a model for an ideal gas. If p denotes its (point) density, the process is characterised 
|13[ [14] [2] by the two requirements that the number of points in a (measurable) set ic 1 
is Poisson distributed with parameter pX(A), where A is Lebesgue measure, and that the 
number of points in sets A\, A2, ■ ■ ■ , A m are independent random variables, for any collection 
of pairwise disjoint, measurable subsets of K . The Poisson process is a model for an ideal 
gas of pointlike particles. 

Let us consider such a process, with density p. Due to ergodicity, almost every realisation of 
it possesses a natural autocorrelation, and the latter can be calculated via the Palm measure 
of the process [2] Thm. 3]. The result reads 

(3) 7 P = P 5 + p 2 X, 

where <5o is the normalised Dirac measure at 0. With 5o = A and A = So, one obtains 

(4) j? = p 2 S + pX, 

which is the diffraction measure. We skip the proof, but mention that, in one dimension, 
the result can easily be derived from the renewal theorem discussed in Section [3] with an 
exponential waiting time distribution [2]. Apart from the trivial point measure at 0, the 
diffraction measure 7^ is absolutely continuous. 
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An interesting modification emerges from a marked Poisson process, where each point of a 
given realisation randomly gets the weight 1 or —1 with equal probability. This leads to the 
following modification of Eqs. ([3]) and dU). 

Theorem 1. Consider a typical realisation of the homogeneous Poisson process of density p 
in M , which is a simple point set A C M . Let uj = ^2 x£j \ W x 5 X be a random Dirac comb 
where (W x ) xeA constitutes an i.i.d. family of random variables that take the values 1 and 
— 1 with equal probability. Then, the corresponding autocorrelation and diffraction measures 
almost surely read 

-f UJ = p5 and ^ = p\. 

Proof. This is the situation of the random weight model of [21 Ex. 7], applied to a stationary 
Poisson process, which is an ergodic and simple point process. The result now follows from 
[21 Thm. 4 and Cor. 1] by a small calculation. □ 

When the density is p = 1, this is one of many examples with diffraction measure A, which 
include the coin tossing sequence on the integer lattice and the Rudin-Shapiro sequence [TT1I5]. 
but also various dynamical systems of algebraic origin [9| such as Ledrappier's shift on Z 2 . 
This provides ample evidence that the inverse problem of structure determination becomes 
significantly more involved in the presence of diffuse scattering. 

One limitation of the Poisson process for applications in physics is the missing uniform 
discreteness. This can be overcome by an additional hard-core condition, as in the classic 
Matern process; see [24] and references therein for a formulation that matches our setting. 
Here, a realisation of a homogeneous Poisson process is randomly marked and then thinned 
out on the basis of a pairwise comparison up to a certain distance. This leads to a modified 
point set that is uniformly discrete (meaning that the minimum distance between any two 
points is a positive number). Despite this modification, autocorrelation and diffraction can 
still be calculated explicitly; see [21 [20] and references therein for more. 

3. Renewal process 

The situation of random point sets is significantly simpler in one dimension, because a large 
class of processes can be characterised constructively as a renewal process. Here, one starts 
from a probability measure p on M + (the positive real line) and considers a machine that 
moves at constant speed along the real line and drops points on the line with a waiting time 
that is distributed according to p. Whenever this happens, the internal clock is reset and the 
process resumes. Let us (for simplicity) assume that both the velocity of the machine and 
the expectation value of p are 1, so that we end up with realisations that are, almost surely, 
point sets in R of density 1 (after we let the starting point of the machine move to — oo). 

Clearly, the process just described defines a stationary process. It can thus be analysed 
by considering all realisations which contain the point 0. Moreover, there is a clear (distri- 
butional) symmetry around this point, so that we can determine the autocorrelation (in the 
sense of ([1])) of almost all realisations from studying what happens to the right of 0. Indeed, 
if we want to know the frequency per unit length of the occurrence of two points at distance 
x (or the corresponding density), we need to sum the contributions that x is the first point 
after 0, the second point, the third, and so on. In other words, we almost surely obtain the 
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autocorrelation 

(5) 7 = 5 Q + v + v 

with v = fj, + fi*fi + fj,*fi*fjL+... and v as denned above, where the proper convergence 
of the sum of iterated convolutions follows from [21 Lemma 4] . Note that the point measure 
at simply reflects that the almost sure density of the resulting point set is 1. Indeed, v is a 
translation bounded positive measure, and satisfies the renewal relations (see [14"1 Ch. XI.9] 
or [21 Prop. 1] for a proof) 

(6) v = fi + fi * v and (1 — fx)v = ju, 

where /2 is a uniformly continuous and bounded function on R. Note that the second equation 
emerges from the first by Fourier transform, but has been rearranged to indicate why the set 
{k | ju(/e) = 1} will become important below. In this setting, the measure 7 of ([5]) is both 
positive and positive definite. 

Based on the structure of the support of the underlying probability measure fA, one can 
now formulate the following result for the diffraction of the renewal process. 

Theorem 2. Let [i be a probability measure on R + with mean 1, and assume that a moment 
of \i of order 1 + e exists for some e > 0. Then, the point sets obtained from the stationary 
renewal process based on /i almost surely have a diffraction measure of the form 

7 = (7) pp + (!"/>) A, 

where h is a locally integrable function on R that is continuous except for at most countably 
many points. It is given by 

_ 2(|£(fc)| 2 -Re(£(fc))) 

m ~ \^m? • 

Moreover, the pure point part is given by 

J 5 , if supp(/i) is not a subset of a lattice, 



(7) , . 

PP l^z/6' tf ^ * s ^ e coarses t lattice that contains supp(;u). 



Proof. The process has a well-defined autocorrelation 7 as outlined above and given in Eq. © . 
Due to the ergodicity of the process, this means that almost every realisation of the process 
is a simple point set with this autocorrelation. Since 7 is positive definite, it is Fourier 
transformable, with 7 being a positive measure on R. 

The point measure at (which is always present) reflects the fact that the resulting point 
set almost surely has density 1; see [H Thm. 1] and its proof for a detailed argument. The 
functional form of h can be calculated from the second renewal relation in ([6]) whenever one 
has fi{k) ^ 1. Its local integrability everywhere is a consequence of the assumed moment 
condition, by an application of |25t Thm. 1.5.4]. 

The distinction via the nature of supp(/i) takes care of the set {k G R | fl{k) = 1}, which 
is either {0} or countable. The result now follows from [21 Lemma 5 and Thm. 1] together 
with Remark 3 of the same paper. □ 
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The renewal process is a versatile method to produce point sets on the line. These include 
random tilings with finitely many intervals (which are Delone sets) as well as the homogeneous 
Poisson process on the line (where \i is the exponential distribution with mean 1); see [21 Sec. 3] 
for explicit examples and applications. 

4. Random matrix ensembles and random point sets on the line 

The global eigenvalue distribution of random orthogonal, unitary or symplectic matrix en- 
sembles is known to asymptotically follow the classic semi-circle law. More precisely, this 
law describes the eigenvalue distribution of the underlying ensembles of symmetric, Hermit- 
ian and symplectic matrices with Gaussian distributed entries. The corresponding random 
matrix ensembles are called GOE, GUE and GSE, with attached /3-parameters 1, 2 and 4, 
respectively. They permit an interpretation as a Coulomb gas, where f3 is the power in the 
central potential; see [1] [21] for general background and |12j for the results that are relevant 
to our point of view. 

For matrices of dimension N, the semi-circle has radius \/2N/ir and area N. Note that, 
in comparison to [2TJ, we have rescaled the density by a factor 1/^/vr here, so that we really 
have a semi-circle, and not a semi-ellipse. To study the local eigenvalue distribution with our 
application in mind, we rescale the central region (between ±1, say) by y / 2N/ir. This leads, 
in the limit as N — > oo, to a new ensemble of point sets on the line that can be interpreted 
as a stationary, ergodic point process of intensity 1; for j3 = 2, see [TJ Ch. 4.2] or [23] and 
references therein for details. Since the process is simple (meaning that, almost surely, no 
point is occupied twice), almost all realisations are point sets of density 1. 

It is possible to calculate the autocorrelation of these processes, on the basis of Dyson's 
correlation functions |12| . Though the latter originally apply to the circular ensembles, they 
have been adapted to the other ensembles by Mehta [21] . For all three ensembles mentioned 
above, this leads to an autocorrelation of the form 

(7) 7 = <5o + (l-/(M))A 



(8) /(r) 



where / is a locally integrable function that depends on j3. Defining s(r) = — ^—*-, one obtains 
(with r > 0) 

' a (r) 2 + S / (r)J r °s(t)dt ) ap = i, 

s(r) 2 , if j3 = 2, 

^s(2r) 2 - 2s'(2r) J Q r s{2t) dt, if /3 = 4. 

The diffraction measure is the Fourier transform of 7, which has also been calculated in 
[12[l21j. Observing 8q = A and A = 8q, the result is always of the form 

(9) 7 = So + (1 - b(k)) X = 5 + h{k) A, 

where b = f. The Radon-Nikodym density h for (3 = 1 reads 

>|(2-log(2|fc| + l)), if|fc|<l, 

2|fc|+l 
2|fc|-l' 

where fegl. The result for (3 = 2 is simpler and reads 
y ' y ' |1, if \k>l, 



(1°) M k ) = S ,,, , 2|fc|+l . f |,| x1 

2 ~l fc l l °gmhr> if|fc|>i, 
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Figure 1. Absolutely continuous part of the autocorrelation (left) and the 
diffraction (right) for the three point set ensembles on the line, with j3 S 
{1, 2, 4}. On the left, the oscillatory behaviour increases with (3. On the right, 
(3 = 2 corresponds to the piecewise linear function with bends at and ±1, 
while (3 = 4 shows a locally integrable singularity at ±1. The latter reflects 
the slowly decaying oscillations on the left. 



while (3 = 4 leads to 
(12) 



h A (k) 



(2-log|l-|fc||), if|fc|<2, 
if |fc| >2. 



Figure [T] illustrates the three cases. To summarise: 

Theorem 3. The eigenvalues of the Dyson random matrix ensembles for (3 £ {1,2,4}, in 
the scaling of the local region around as used above, almost surely give rise to point sets of 
density 1, with autocorrelation and diffraction measures as specified in Eqs. ([7]) and ([9]). □ 

Note that /14 is smooth at k = ±2, but has integrable singularities at k = dbl. The latter 
are a consequence of the stronger oscillatory behaviour of the function f± at integer values, as 
was already noticed in [12]. When extrapolating to other values of (3 (in particular to j3 > 4), 
this is the onset of another Bragg peak. 

It is well-known that the circular random matrix ensembles (COE, CUE, CSE) asymptot- 
ically give rise to the same local correlations [121 EH j an d hence to the same autocorrelation 
and diffraction (after appropriate rescaling). 



5. Random matrix ensembles and random point sets in the plane 

The above examples were derived from matrix ensembles with real eigenvalues, and thus 
lead to point processes in R. There is also one ensemble, due to Ginibre [15] (see also [H]), of 
general complex matrices with Gaussian distributed entries that will give rise to a stationary 
point process in R 2 . Again, this emerges (by proper rescaling) from the eigenvalues (now seen 
as elements of the plane), which approach uniform distribution in a circle of radius ^/N/it 
(and hence area N) as N — > 00. 
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Figure 2. Radial dependence of the absolutely continuous part both of the 
autocorrelation and the diffraction measure for the planar point set ensemble, 
as derived from Ginibre's matrix ensemble. 



As before, the system can be interpreted as a Coulomb gas, with a potential parameter 
/3 = 2. Other matrix ensembles permit this interpretation, too, but do not seem to correspond 
to interesting stationary processes, wherefore we stick to Ginibre's example here. 

Following the original approach of [15], the limit N — >• oo leads to a stationary and ergodic, 
simple point process of intensity 1, so that almost every realisation is a point set in the plane 
of density 1. Using complex variables Zi G C ~ M. 2 , the 2-point correlation function is of 
determinantal form, 



(13) 



p(zi,z 2 



e -K{\z 1 \ 2 + \z 2 \ 2 ) 



= 7r l 2 ll 



=*"l z 2 



'1 



-■k\z x -z 2 \' 



)• 



see [15] or [21] for a derivation. Note that, despite using complex coordinates here, the 
expression is calculated relative to the volume element of real coordinates (hence relative to 
Lebesgue measure, as in [21]). The result is translation invariant and only depends on the 
distance r between the two points. 

As a consequence, the autocorrelation of a realisation almost surely reads 



(14) 



7 = 5o + (l-e-^)A, 



which is radially symmetric, with r as above. By a standard calculation, the Fourier transform 
of 7 results in 



(15) 



7 



<*o + (l 



-l fc l 2 )A, 



so that we obtain a self-dual pair of measures under Fourier transform (as in the Poisson 
process of density 1). The radial dependence is illustrated in Figure EJ 

Theorem 4. The Ginibre complex random matrix ensemble, in the scaling used above, almost 
surely results in point sets of density 1, with autocorrelation (|14p and diffraction (|15p . □ 
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6. Summary and Outlook 

In this short communication, we have discussed several explicit examples of stochastic point 
sets with explicitly computable autocorrelation and diffraction measures. The viewpoint of 
point process theory provides a universal platform to do so, though our examples above also 
admit a direct approach. It would be interesting to extend this to a family of processes, 
with /3 as parameter in the spirit of |1 H I27 | IT5]. which then interpolates between the Poisson 
process of density 1 (/3 = 0) and the integer lattice (which is approached as )3 — > oo). 

An interesting question concerns the connection with dynamical systems theory, in partic- 
ular the general relation between diffraction and dynamical spectra. Recent progress suggests 
that such a connection might indeed exist, although it will certainly be more involved than 
in the pure point diffractive case. 

One fundamental shortcoming so far is the lack of understanding and explicit examples 
for randomness with interaction. A first step in that direction needs the inclusion of Gibbs 
measures for equilibrium states, though it is not clear at the moment to what extent one can 
derive explicit examples (such as the classic and well-known Ising lattice gas). 
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